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Abstract 

The formation energy of a sohton was derived in the TLM model to be 2Ao/vr which 
is independent of the masses of the polymer atoms. We estimate effects of the electron- 
phonon coupling on the soliton formation energy, particularly, through modifications of the 
harmonic oscillations of (CH) groups. A model based on the amplitude mode formalism is 
used and the electron-phonon coupling is taken into account to make the effective potential 
of lattice vibration. Among dynamical degrees of freedom in the lattice vibration, a trans- 
lational mode of the soliton is separately examined, since it corresponds to a zero-energy 
mode. We get a difference in the formation energy between the two isotopes of about 
0.03eV. Soliton trapping around defects or impurities changes the formation energy and 
enlarges the isotope difference. 

§1 Introduction 

The trans-polyacetylene is an ideal material which has excitations such as solitons, polaron, 
bipolaron, and so on. They can be described very well by the SSH model which was 
proposed by Su et al% This model contains the electron nhonon interaction in addition 
to electron hopping and lattice vibrations. Takayama et aZBP made the TLM model which 
is a continuum version of the SSH model. It makes the analytical study of polyacetylene 
possible. 

A number of experimental studies have supported the presence of a kink excitation in 
t-polyacetylene {t-{CYi)x) ■'^^'^'^ Suzuki et al. have reported an optical absorption measure- 
ment in lightly doped polyacetylene .0) They have found that the absorption related to the 
electronic mid-gap state is enhanced by doping. This result has been interpreted as tran- 
sitions between the mid-gap levels of solitons, induced by the doping, and the conduction 
band. 
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In addition to the absorption due to the mid-gap scattering transition, sevral absorp- 
tions have been observed near 1370cm~^ and 900cm~^ in the experiments on the infra-red 
absorption for a hghtly doned (CH)^;.^ Mele and Rice have suggested that these peaks 
are due to lattice vibration.S^ Horovitz has explained these peaks as the lattice vibration 
around a soliton added to a dimerized state, proposing a model with several phonon bands. 
This model is called "the amplitude mode formalism" 

On the other hand, the experimental studies of the photo-induced absorption spectra of 
undoped t-(CH)^ have also shown additional infra-red active modes. The positions are 
listed in Table I. These infra-red active peaks have been considered to be due to internal 
vibrational modes created by the presence of the soliton. The intensities of these new modes 
are much weaker than those of the first three modes. 

Ito et al. have studied the phonon structures in the presence of a soliton in the frame- 
work of the continuum theory for trans-polyacetylene (TLM model). They have analyzed 
small lattice oscillations by using time-dependent adiabatic mean-field equations resulting 
from a selfconsistent treatment of the electronic field and the phonon field.B* They have 
found three localized modes: the Goldstone mode (or the translational mode), the ampli- 
tude oscillation mode (the shape mode) and a third mode. Ito et al. have also found that 
two of the three localized modes, i.e. the Goldstone mode (G.M.) and the third localized 
mode are IR active, whereas the amplitude oscillation mode is inactive.ll3) 

Mele and Hicks have applied Horovitz's amplitude mode formalism to investigate the 
optical properties of the localized phonon modesM They have shown that there are three 
absorption peaks due to the translational mode in the presence of a single soliton. Terai, 
Ono, and Wada have applied the Mele-Hicks theory to calculate the phonon structures, 
taking account of the translation and the third localized modes, where they have used three 
phonon bands.ll3 (For each lattice site, there are several degrees of freedom such as bond 
vibration between C and H and that between nearest neighboring C atoms. Within the 
polymer plane, there are 6 degrees of freedom.ll3) The number of IR-active modes e.g. gerade 
modes is, explicitly, 4.) Heeger et al. performed an optical-absorption experiment with high 
resolution, observing two addtional IR peaksJ13' Terai et al. have estimated relative optical 
absorption intensities and obtained a good agreement with the experiment JIM^ 

In this paper, we would like to study an isotope effect on soliton formation energy in 
t-polyacetylene. attice oscillations .Ej^ 

of the unharmonicity of the effective potential gives a correction to the phonon frequency 
and the electron-phonon coupling constant. We will show that the correction can give the 
difference in the formation energy as large as the observed value. 

§2 The effective potential of the one-dimensional electron-phonon system 

In quasi-one-dimensional electron phonon system bearing the Peierls transition, there are 
two types of phonon modes. One is the amplitude mode and the other is the phase mode 
(the CDW sliding mode). In the incommensurate charge density wave, the ion displacement 
at the mth site is 

Um = uq cos(qo fna + 5) 
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where uq is the displacement amplitude, qo is the wave number which is equal to Ikp- The 
phase variable 6 describes the motion of the charge density wave. In the half-filling case 
such as trans-polyacetylene, qq = ir/a, so that Um = {—)^uqcos 5. Thus the dynamics of S 
in t-polyacetylene is written by using only the amplitude of the charge density wave. 

Let us consider a one-dimensional electron system including a electron-phonon coupling 
with a optical phonon field A = quq where g is the electron-phonon coupling constant, the 
Hamiltonian is given by 

+ J2 I (^x^\{x)[-WF'J^d^ + <yi^{x)]'4>s{x), (2.1) 

s •' 

where ujq is the bare phonon frequency. The electron field V's is introduced as a two 
component spinol field representing left and right-going waves at the Fermi level. A is 
a dimensionless electron-phonon coupling constant, the Hamiltonian ( |2.1| ) is called the 
Takayama-Lin-Liu-Maki (TLM) model. 0) 

The coupling between electron and phonon fields induces unharmonicity of the lattice 
vibration and, due to this non-linearity, the Hamiltonian ( ^.l| ) bears a bond order wave 
state (BOW state). 

To investigate the unharmonicity, we make a eff'ective potential for the phonon field 
using loop expansion of a partition function. 

The partition function of (2.1) is written as follows: 




Z = Zo ( Trexv - / ^/dr ) . (2.2) 



where (Tr) being a thermal average with respect to the unperturbed Hamiltonian. Hi is 
the perturbation term given by the electron-phonon coupling '\\)\a\l\'\\)s. Zq is the partition 
function of the unperturbed Hamiltonian in which the electron-phonon coupling is excluded. 
Zq is written in the path integration form, 



Zo = j{D^D^l,lD^l,s)e^v[- j '^r{L\ + L%)^ 



dr 



4 = ldx4ix)^tM + Y.I<^^^l(^)i-^^Fa3dMx), (2.3) 



The thermal average in ( |2.2[) is also written by the path integration form. The path integra- 
tion over tl^s is approximated by a diagramatic calculation. Consider the summation over 
connected loop diagrams which are shown in Fig.l. The solid line is associated with 
the electron Green function. The solid circle corresponds to the vertex with ai coupling. 
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In the momentum space, the electron Green function is given by 



dxdr (Tri)s{x,T)4)l,{x' ,0)) e 



'ikx—iujT 



iwl + vpka^ 

WW 



(2.4) 



where s denotes up or down of the spin orientation. Representing the sum of the diagrams 
in Fig.l as U, the partition function is approximately given by 



Z 



/(DA)exp< 


-JdTLl + Jdrdxu'^ 


/(Z)A)exp< 




\ 



(2.5) 



Let us estimate U. The lowest order term with respect to the momentum transfer q between 
the electron, ^g, and the lattice, A{x), is given by the summation written as follows: 



1^ n..<^.iji 



dk- 



— [A{0)aiG^^\iLO,k)Y 



4:1171 



^ V / dfcln (l 



A(0)^ 



(2.6) 



where the summation over uj is taken over the Matsubara frequency lo = T{2m + l)7r, m 
being an integer. A(0) is the Fourier transform of A{x) at = 0. Taking zero temperature 
limit and using the relation 



1 v-^ , fiuj — \a\ 
oj.odd 



1 
47ri 



de tanh — In — 
2 e 



1^1 



1 



(\a\-\b\), (T-O). 



we can rewrite 02. 6f) as 



TTVp 



-A\0) In 



m 

2A 



(2.7) 
(2.8) 



where A is the electron band cut off given by A = vp/a- 

The higher order terms with respect to q is given by diagrams shown in Fig. 2. The 
summation over these diagrams gives 



LiJ,S 



ILO + Vpik + g)(73 + A(Ji 

(ia;)2 -4(A: + g)2 - A2 



A(Ji 



iuj + vpka^, + A(7i 
(iw)2 - {vpk)'^ - A2 



Acji 



iw + vpkas + A(Ji 
(ia;)2 - [vFkf - A2 



-A 



TTVp 



v/4A2(0) + {vpqy^ f vpq+ ^ 4A:^{{}) + [vpq)^' 



vpq 



2A(0) 



A. 



(2.9) 
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where we subtracted from the first term in the right hand side of ( |2.9|) the q independent 
term (the second term). In the real space, this term would be expressed as 



1 



■A(x) 



7 



ln(7 + Jl + 72)-l A(x) 



where 



7 



VF 



d 



(2.10) 



(2.11) 



2A(0) dx' 

A(0) is considered to be absolute of the mean value of A(x). Since U is expressed by the 
sum of (p. 81) and (|2.9D, we can study the lattice dynamics of the electron-phonon system 



(2.1) by using the total effective potential which is given as follows : 



1 



dx 



X 



+ 2A(x) 



7 



ln(7 + Jl + 72) - 1 A(x) + 2A(x)^ In 



|A(x) 



A{xf 



(2.12) 

In the perfectly dimerized state, and in the adiabatic limit, the spatial variance of A{x) is 
equal to zero. The effective potential have double minimum points at A = =bAo where Aq 
is determined by 



Ao = 2A exp 



1 

2A 



(2.13) 



The potential is written by the dashed line in Fig. 3. Thus, we consider that |A(0)|, which 
was introduced in ( |2.10| ) as the mean value of A(x), is given by Aq. Therefore, the effective 
Hamiltonian is given by 



H 



2ttvfXu>q 
1 



dx 



+ 



2ttvf 



dx 



A{xf+u^lA{x)\ 
2A{x)D^A{x) + 2A{xf 



In- 



\A{x) 

1a" 



(2.14) 



where A{x) is the time-derivative of A(x). and 7 are introduced as follows 



7 



ln(7 + Vl + 72)-l, 
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(2.15) 
(2.16) 



' 2 dx' 

where ^0 = vf/Aq. As is well-known, the electron-phonon system (|2.l| ) bears the non- 
linear excitation. A soliton is the solution of ( |2.1| ) which seperates A(x) = Ag phase and 
A(x) = — Aq phase. The effective Hamiltonian ( p.l4| ) also has a soliton solution. It's spatial 
form is shown in Fig. 4. We obtained this form by using the variational method. (Appendix 
A.) The form of A(a;) for the soliton is nearly equal to the well-known form. 



A<j(a;) = Aotanhx/^, 



(2.17) 
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pee 
|) In 



which is the solution of the Hamiltonian Our result gives a value of the soliton 

width of about ^ = 1.04^o with = vf/^o- 

The soliton excitation energy which is centered around the soliton is mainly kept in the 
electron system. In a soliton case, the energy level of the electronic states in the va^ 
band are shifted by the lattice distortion and the mid-gap localized state is induced, 
our model, the property of the electron system is approximated by the effective potential 
which is constructed on the assumption of the continuum electronic dispersion. Since the 
localized electronic state is not included, our model, ( p. 14 ), does not include the exact 
change of the electronic energy level. Using (|2.17| ), and using the fact that As{x) is the 
solution of ( [A.ip , the soliton formation energy is derived to be Aq/vl^ (Appendix A). Our 
result gives the different value with the result given by the TLM mode. 10^ On the other hand, 
the Hamiltonian ( p. 14 ) gives good coincidence with the TLM model, (2.1), in analyzing the 
lattice dynamics of polymer chain. As we will show later, the dispersion relation of phonon 
and the localized ocillation mode around the soliton obtained in the TLM model are well 
reproduced in our model. Then, the Hamiltonian ( |2.14| ) would be suitable for analyzing 
the renormalization of the unharmonic lattice oscillation and for making the correction of 
the parameters A, wg, etc. 

The eigen mode gi of the lattice fluctuation rj around a static soliton is given by the 
Schrodinger equation 



2Xiu; 



-91 



ln(7+ Jl + 72) - 1 



+ 1 + In 



|A,(x) 
An 



(2.18) 



where As(x) is the static soliton solution of ( |2.14| ), i.e. ( p. 17 ). In the perfectly dimerized 
state, the log term of the potential in ( p.l8| ) vanishes and the continuum spectral of phonon 
frequency $7^ is given by 



sinh s, 



(2.19) 



where s = vpq/ {21S.q), q being the wave number of the phonon. Then, the extended modes 
of lattice fluctuation around the soliton have the spectrum like ( 2.19| ). This dispersion 
relation was also derived by Nakahara and Maki using the random phase approximation 
(RPA) method for the electron polarization.c3) 

In the Hamiltonian (p.l|), there are localized modes which have non-zero amplitude 



around the soliton center. They are the Goldstone mode, amplitude oscillation mode and 
the third localized mode of soliton. These localized modes have been obtained by several 
groups using the RPA method. 0)®- "ES We show in Appendix B that ( p. 18] ) has the zero- 
energy mode (the Goldstone mode) and that its form is given by Substituting the 
form of the soliton which is shown in Fig. 4 into ( 2.1^ ), the frequencies of the other localized 
modes are numerically obtained. (ApnendbL B) Their frequencies are shown in Table II. 
Our result coincides with other works .t3)E3)EfP 

In this paper, we discuss the quantum correction to the soliton formation energy using 
(|2.14| ) and soliton solution (|2.17D . 
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§3 Anharmonicity of the lattice vibration 



Because the potential of the lattice field A has the logarithmic term, the lattice fluctuation 
r] around A = zbAg has an anharmonic potential. In ( |2.1| ), the parameters ojq, A and vp 
were determined so as to obtain good agreements with experimental values of the electron 
band gap and the phonon dispersion considering the lattice fluctuation as the harmonic 
lattice oscillation. In the dispersion relation ( ^.191) and the selfconsistency equation for the 
electron band gap, ( p. 13 ), these parameters were used. However, since anharmonicity of 




the lattice oscillation modifies the mean position of (CH) groups and the energy spectral 
of phonon dispersion, we should make corrections for the parameters A, ujq and vp. In the 
perfectly dimerized case, expanding A{x) around A = Aq, the Hamiltonian becomes 

1 f ( fp' 2A 

where A(x) = Ag + r/. Suppose that the higher order terms in r]/ /S.q are sufficiently small in 
the quasi-one dimensional system, we neglected terms over 5th order with respect to r//Ao. 
We compare the potential of ( ^ ) and (g) graphically in Fig.3. The potential of ( pO| ) 
is drawn by the solid line. The coincidence is fairly well when |r//Ao| < 0.5. For the realistic 
material, t-polyacetylene, the parameters are A ~ 0.2, loq ~ 0.3ey, and Aq ~ 0.67eV, then 

^ = (T,7?(x, r = 0+)7?(x, 0)) = — 2 / d{vpq)^ ~ 0.4, (3.21) 



Ag ^ 4Ao 



where ujq is given by ( 2.19 ). Considering that over the 5th order contribution from rj/lS., 







is negligible, we investigate the effect of lattice fiuctuation described by ( 3.20| ). The bare 



phonon Hamiltonian which is written by the first three terms in the right hand side of 
( 3.20| ) gives the bare Green function of ry. 



Dr,{uj,k) = - {Trr]{x,T)r]{x' e 



-ikx ^—iujT 



irvpXujQ 



- r ,2 02- (3.22) 

The forth and the fifth terms in the right hand side of ( |3.20| ) give 3 point and 4 point 
vertices, respectively. They are associated with the diagrams in Fig. 5. 

§4 Correction to the soliton formation energy 

The 3 point and the 4 point vertices give the correction to spectral of the harmonic oscil- 
lation. We will make one-loop approximation and estimate the correction to the soliton 
formation energy. In the one loop approximation, the 3 point vertex gives the effective 
contribution which is shown by diagrams in Fig. 6a). 

/ dx?7 • - • — (r?\ = / dx — ^5m ri. (4.23) 

J 2 TTVfAq \ I J TTVp 
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In the same way, the 4 point vertex gives the effective contribution which is shown in 
Fig.6b). 



dxrj 



1 



2'nvp 



-5m. 



(4.24) 



Here, we introduced the notation 5m = (rf) /^o- Totahy, the contribution becomes 

5m 2\ 



2'kvf 
5m 



j dx(-2?7Ao + r]^ 



2-KVF 

5m 
2ttvf 



dx 



dxA^. 



i-l)A(x)2 + 2A(x)2ln^ 



(4.25) 



Since r/ = A(x) — Aq is a small quantity, and since the functional form of the correction 
term should be similar to the effective potential in ( |2.18D , we approximated the first line of 
(|4.25| ) by the second line. 

Because the parameters A, ujq, and in ( ^Tll ) is selected to give the good coincidence 
with the realistic material in the mean-field theory, the Hamiltonian should be reconstructed 
to cancel the contribution ( 4.25| ). As the corrected Hamiltonian, we obtain 



1 



2itvf\oj'q 
1 



2'nvF 



dx 



2A{x)D^A{x) + 2A{xf In - A^{x) 



2A 



+- 



+- 



5m 

2'KVF 

5m 

2'KVF 



dx 



All. 



1 A(x)^ + 2A(x)2ln 



|A(x) 
2A 



A2(x) 



(4.26) 



where 7 = — i^^^- The third and the forth line in ( 4.26| ) is called as the "counter term" 
which is added to cancel the contribution ( [4.25| ). In the derivation of ff^, we assumed that 
the fluctuation ry/Ag is smaller than 1, and terms over the 4th order with respect to tj/Aq 
is negligible. In the quasi one-dimensional electron-phonon system, like polyacetylene, this 
condition would be satisfied. In the purely one-dimensional system, the higher order terms 
with respect to ??/Ao is not negligible. Then, we can not renormalize the anharmonicity 
into the form of ( |4.26D . We do not consider the purely one-dimensional case. 

Next, we renormalize the unharmonic contributions which come from the counter terms 
in ( [4.26 ). Repeating these procedures, we obtain the renormalized Hamiltonian which 
represents the actual situation of the quasi one-dimensional electron-phonon system at 
the half- filling. By the redefinition of the parameters A, vf and ojq, the renormalized 
Hamiltonian is written in the form of ( p. 14 ). Comparing ( |2.14| ) and ( [4.26] ), we obtain the 
renormalized Hamiltonian 

1 



Hr 



2kvfrXrujq 



R 



dx 



ARixf + iOQji^A 



R 



(xY 
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1 



2TTVFR 



dx [2dRAR{x)D^AR{x) 



(4.27) 



(1 - 5m)-^/^{l + Aln(l - 6m)y/^ujQ 
= {1 + Xln{l - 6m))-^ X 
Ar = {l-5m)-^/^A 
dR = (1 - 5m) 

VFR = VF 

Ar = VFR/a (4.28) 



where a is the lattice constant. We determined ( 4.2g| ) on the condition that vfr and Ar 
are invariant under the renormahzation since they are parameters of the electron system. 
Because 6m > 0, the corrected parameters ^^Qpj and \r are larger than the bare values. 

From ~ {^—6m)~^ , we can see that, in the actual material, the effective potential 

is (1 — 5m)~^ times deeper than the uncorrected one. 

On the other hand, since the effective potential becomes deeper, the differential equation 
for the soliton is affected by the renormahzation. We calculated the soliton width £^r for 
the various values of 1/dR = (1 — 5m)~^ using the variational method. (Appendix A) We 
plot £,o/£,R. in Fig. 7 as the function of 1/dR. The soliton width becomes narrower as the 
effective potential becomes deeper. It behaves like 

^ = 1.04dR°-^2 ^ 104(1 - 5m)°■^^ (4.29) 

in the range 0.2 < dR^^ < 2. 

In the static limit, A = 0, using ( [A.l| ) , the soliton formation energy is given by 



— / dx{Aol - Ar{x) 
Vfr J 



2'kvfr 



^ (4.30) 

where, from the first to the second line, we used the soliton form A(x) = Aq tanh(x/^/j). 
In addition to ( |4.3C| ), zero-point lattice oscillation gives a contribution to the total energy, 
which is explained by the path integration.0) The presence of the soliton changes the 
phonon frequency and this contributes to the soliton formation energy. Then, the soliton 
formation energy is, totally, given by 



n ^ soliton f, 



dimer^(l-M°-^^ vr"" ^^"^^^ 
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§5 The effect of impurity 



In the reaUstic material, there are lattice defects or impurities. The potential due to the 
impurities traps the soliton and therefore the formation energy is affected by the trapping 
energy. Suppose spin or charge of the soliton interacts with impurity and the impurity 
potential is given by the delta function —V5{x). If the soliton is rigid, the trapping energy 
will be 

no 



V / dx5{x)n{x — X) 



-V- 



cosh^ X/e' 



(5.32) 



where X is the position of soliton center, V the magnitude of the impurity potential, ^ is 
the soliton width, and n{x) is the charge or spin density of the soliton. We approximated 
n{x) as rinj cosh^(x/^) because the electron density in the mid-gap state of the soliton has 
this form.cP Since we are interested in the dynamics of the soliton center and the effect of 
soliton trapping. We use the effective Hamiltonian. 



tr 



2M 



U 



cosh^X/C 



(5.33) 



with U = \ VnQ\. Here M is the soliton mass which is given by 



M 



2'nvF 



dx 



1 dA(x) 
K'n dx 



(5.34) 



and P is the momentum conjugate to the position of the soliton center X. ( p. 33 ) gives the 
trapping energy spectrum given as follows: 



En 



8Me 



1 + 8Mf7^2 _ (2n + 1) 



(5.35) 



where n being the non-negative integer which satisfies < n < {^Jl + ^MUS^ — l)/2. The 
experimental value of the soliton formation energy would correspnd to the sum of "pure" 
formation energy of the soliton Eg and the trapping energy Eq. The formation energy of 
the soliton is lowered by 



En 



1 



8Me 



1 + 8MU^^ - 1 



(5.36) 



In the weak trapping limit, MU^^ « 1, Eq becomes -2MU^^^. When MU^^ » 1, Ei 
becomes -U + ^U/{2M^^). 



§6 Isotope effect 

The realistic t-polyacetylene has several normal modes of lattice oscillations. Horovitz has 
proposed the model which include these normal modes introducing plural components of 
phonon fields. 
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The multi-mode Hamiltonian is given 

3 . 

y - 



H 



dx[Ai{x) + LviAiix 



(6.37) 



where 3 types of phonon fields are taken into account, ujq being the renormalized phonon 
frequency 



^0 = 2Aw^. 



A(x) is defined by 



.38) 



where Aq is the coupling constant of the ath phonon with the bare frequency oja- There is 
a relation 

A = 5^Aa. (6.39) 

a 

Using ( |6.37| ), the effective potential is obtained as in the same way as shown in the section 
2. Introducing a phonon field of the ath normal mode by rja, A is expanded around A(a;) 
as follows: 



A{x) + r/, 

Aa UJa 
A iOo 



■Ha- 



A = 

a 

Substituting ( |6.40 ) into the effective Hamiltonian ( |2.14| ), we obtain 



(6.40) 



+— y / dx 

a,l3 



—Va-L>yV/3 



X W, 







+ 1 + In ■ 



A(x)| 1 \ y^XaX/s UJaUJfB 



Ac 



2A 



A 



VaVf3 



(6.41) 



where is given by ( |2.15 ). The eigen mode of the Hamiltonian is constructed as follows. 
From ( |6.41 ), the Schrodinger equation of an eigen mode (fa is given by 



\Aix) 



+2 ^ J XaXf3LJaUJ/3 1 + In 



An 



2a' 



(6.42) 
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Since we are considering the one-dimensional case, the eigen-mode function of ( |6.42| ) can 
be expressed by the super position of 

^a = y2^a/9£, (6.43) 



where gc is the eigen function of ( p.l8 ). We write the frequency of ge by O^. Substituting 



( 6.43| ) into ( |6.42| ), and using eg. ( p. 181 ), Ca^pe is given by the equation. 



This equation was obtained in a different way using the RPA method by Terai et al^^ 
From i ^M ), the expansion coefficient Ca,e and the frequency is obtained. J7 is given by 
the solution of the equation 



2 



-Dm, 



a 



There are three solutions for each £. Let us call them as Qp^i {(3 = 1,2,3). Using ( |6.44 ) 



and (6.45), Ca/ is given by 



Using (|1|), (|6l5|) and (|6l6|) , the ph onon frequency and the eigen-mode function are 
obtained. 

Substituting the values in Table II into ( |6.45| ), the localized modes around a soliton 
are obtained. They are summarized in Table I for isotopes of (CH)^. and (CD)^,. The 
experimental results are also shown. In the calculation of the phonon frequency, we used 
the parameters and Wq, shown in Table III. They were determined by Vardeny et al. 
from the Raman scattering experiment. The phonon frequencies in (CH)^ are relatively 
larger than those in (CD)^. 

Let us compare the formation energy of the soliton between the two isotopes of t- 
polyacetylene. As we have discussed in §3, the renormalized Hamiltonian is constructed 
by adding the counter term, and the soliton formation energy should be corrected by the 
fluctuation Sm = {rf') /Aq. Using the dispersion relation Vt^^k which is obtained from ( |6.45 ) 



for the dimerized state, and using the phonon Green function which is derived from ( |6.41 ), 
5m is estimated as follows : 



5m = (Tt r = 0+)7?a(a;, 0)) = ^ / d(uFg)^^ 

/ 0.4536 for (CH), 



0.4002 for (CD), 
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In the calculation in ( |6.47| ), we used Qa,q which is given by ( |2.19| ). In Appendix B, we have 
solved ( p.l8 ) for the soliton approximately. Since we have made an expansion of up 



to 6th order with respect to 7, we cannot exactly reproduce the dispersion relation ( p. 19 ). 



Then, we will use the phonon frequencies given by the RPA method which is explained in 
ref.jl^. Substituting the phonon frequencies and ( 6.47| ) into (4.31), the correction to the 
soliton formation energy is totally given by 

/ 0.0323eV, for (CH), 
- { 0.0033eV, for (CD).. ^^'^^^ 

On the assumption that the difference in the electron band gap Aq between the two isotopes 
is very small, the corrected parameters are obtained by ( |4.28 ) where the parameters A = 0.2 



and Aq = 0.656eV are used for both isotopes. The corrected coupling constant A/j for(CII)a; 
becomes larger than that for (CD)^;. Then, the soliton formation energy in (CH). becomes 
about 0.029eV larger than that in (CD)^;. 

The isotope difference of the solito n fo rniation energy is observed by the photo-induced 
absorption experiment by Schaffer et aZ.lilPta) They concluded that, from the observed value, 
the difference is about 0.03eV ~ 0.04eV. This is of the same order with our theoretical es- 
timation. However, the experimental result is a bit larger than ours by about O.OleV. This 
difference would come from the trapping by impurities. In the optical absorption experi- 
ment, the frequency of Ti mode, which corresponds to the excitation of soliton translational 
motion in (CH). is larger than in (CD)^;. This difference is considered to come from the 
soliton mass Mg. In the multi-component theory (the amplitude mode formalism), the 
soliton mass is given by c3i 

ttvfX'^XujJJ \ dx J a ^\K^J J dx 

(6.49) 

where Ma is the effective mass of the ath lattice oscillation mode and Ka is the effective 
spring constant of the ath mode. They are introduced by the following relation. 



Ao = guo, 

A : 



Since the a bonds between atoms get small influences from the hydrogen mass, the 
isotope difference of would be small. On the other hand, the effective mass of lattice 
oscillation in (CH). is smaller than that in (CD)^;. Therefore, from ( 6.49| ), the soliton mass 



in (Cll)a; is smaller than in {CI})x- Thus, in sufficiently deep impurity potential, from 
(|5.36| ), the trapping energy in (CD)^; is lower than (CH)., and the excitation energy from 
the lowest bound state to higher excited state in (CII)^; is larger than that in (CD)^^. In §5 
we saw that the soliton formation energy is affected by the soliton trapping energy. Suppose 
that the impurity potential is sufficiently deep and that the potential is well-approximated 



13 



by the harmonic potential at the bottom. Considering the soHton trapping energy as the 
energy of the zero-point oscillation around the impurity potential, the isotope difference of 
the trapping energy is about (0.062 — 0.0507)/2 = 0.006eV, where the experimental data 
of Ti mode which is shown in Table I is used. Therefore, the soliton formation energy in 
(CH)^, is larger than that in (CD)^, by 0.0029 + 0.006 = 0.035eV. This theoretical estimate 
agrees well with the observed value. 



§7 Summary 

We investigated the effect of lattice fluctuation on a soliton and compared our result with 
experimental result observed by the photo-induced optical absorption experiment. First, 
we renormalized the electron-phonon interaction into the effective potential of lattice dis- 
tortion field. Then, the effective Hamiltonian is solved to give the soliton solution. We 
obtained the dispersion relation of the lattice fluctuation around the soliton using the effec- 
tive potential and reproduced the theoretical result given by using the RPA method. There 
is a good agreement between our result and experimental result about the level position of 
the localized mode around the soliton. 

Due to the anharmonicity of lattice vibration, the phonon frequency and the electron 
band gap in the actual potential is affected from the mean-field value. We proposed the 



renormalized Hamiltonian ( 4.27 ) and estimated the correction to the soliton formation 
energy. The magnitude of the correction is estimated by the one-loop renormalization 
5m = (rf') /Aq where rj is the lattice fluctuation field and Aq is the observed value of the 
electron band gap. We found that the actual depth of the double-well potential is wider 
than the depth expected by the mean-field theory. The correction to the formation energy 
due to the electron-phonon coupling is given by 

^E, = Y.\^^ , -Y.\^^, +((l-M-°-«'-l)-- (7.51) 

^ soliton ^ ^nmoT- V / TT 



dimer ^ ' 



On the assumption that there are differences in the electron-phonon coupling constant and 
the phonon frequency between the two isotopes, we found that the formation energy in 
(CH)^; is larger than that in (CD)^,. by about 0.029eV. In the realistic material, the soliton 
formation energy would be affected by the soliton trapping by the impurities. We have 
shown that the isotope difference of the soliton formation energy is widened by the trapping 
energy. In the experimental result, the difference between the two isotopes is about O.OSeV 
~ 0.04eV. The difference from our result is of the same order with the trapping energy. This 
suggests the relevance of further theoretical and experimental investigation of the soliton 
dynamics. 
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A. Soliton solution 



The Hamiltonian ( |2.14| ) has a kink solution (soHton) connecting the two minimum points 
A = ±Ao. Let us calculate a static soliton using ( p. 14 ). The Lagrange equation of ( p.l4 ) 
in the static limit is given by 



\ |A| A 

1 A + 2AlnU + _ 



0, 



(A.l) 



where 7 = — i^^- ■^o is given by ^0 = vp/ ^q- Using ( 2.13| ) and introducing the notation 
u = A/Aq, (|A.1J ) is rewritten as follows: 



6eo^ = 2nlnH. 



(A.2) 



Here we expanded the left hand side of ( A.l ) up to the second order with respect to 7. The 
soliton solution is obtained under the condition 



u{x) 
n(0) 



±1, 
0. 



at(x 



±00) 



Multiplying ^ to (| 



and making an integration over x, ([A.2|) becomes 



du 
dx 



^su ( , 7T-, dx = / dx 

Ve Jo ^l + u^{2ln\u\ - 1) Jo 



(A.3) 



(A.4) 



From ( [A.4| ), the soliton is given by 

A(x) = 

giu) = 



du 

y^l + u^{2ln\u\ - 1)' 



(A.5) 



The dotted line in Fig. 4 shows A(x) given by (|A.5D . It is nearly equal to Aq tanh(2.2x)/^o 



(A.5) is narrower than tanhx/^o- However, considering the higher order derivatives in 
( A.l ) and using the variational method, it is shown that the soliton width becomes wider 
and its form is similar to the form tanh(x/^o) • Since the function ( [A.5| ) would be well 
approximated by sgn(z)Ao(l — exp(— 3.88|z|))(l — 0.497exp(— 3.88|z|)), where z = x/^^q, we 
use the trial function 



As{x) = sgn(z)Ao(l - exp(-a|z|))(l - 0.497exp(-a|z|)), 



(A.6) 



where a is a real number. Substituting ( |A.6| ) into ( |2.14| ) and taking the static limit, the 
mean value of ( |2.14 ) is larger than the dimerized state by 



^ -Re 

7Ta \ 



1.249 



— 1 a \ / vl — ^ 

sin"^ - - 1 + 0.1245 sin"^ a 

a 2 \ a 



+ 0.9452 



(A.7) 
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where we have omitted the imaginary part of the mean value. We used the following 
relation. 

A(x)^A(x)dz = -2a2"-^ [0.6305 + 0.12351(2)2"-!]. (A.8) 

(A. 7) has a minimum point at a = 1.70. Then, ( |A.6D is 2.28 times wider than the function 
( A. 5 ). Solid line in Fig. 4 shows ( A. 6 ) with a = 1.70. Its form is nearly equal to the form 
of 

As{x) = Aotanh(3;/0, (A.9) 



with 



e = 1-04^0. 



(A.IO) 



Using that A(x) is the solution of (A.l), the soliton formation energy is given by 

^l-Jdx{Ao'-A',{x))- (A.ll) 



Since the soliton form ( A.9 ) has the same form as the solution of the original Hamiltonian 
(|2.1|) , it is adequate to use As(x) = Aq tanh(3;/(^o) for the estimation of ( A.ll|) . Therefore, 
( A.ll| ) becomes Aq/vt. 

The correction due to the renormalization of the unharmonic potential gives the effective 
Hamiltonian (4.26). When dji is varied, the soliton width ^r, which is obtained by using 
the trial function ( A. 6 ), changes. We plotted for various values of dR~^ in Fig.7. 



B. Lattice fluctuation around a soliton 

Differentiating (SH) with respect to x/(,o = kx, we obtain the following equation. 

dA,(x) , A , , |A,(x)|\ dA,(x) 



7 



ln(7 + Jl + 72)-l 



d{K.x) 



+ 1 + In ■ 



I An 



d{Kx) 



0. 



(B.l) 



Comparing (|B.l| ) with ( 2.1§| ), we find that 
Goldstone mode. 



is an eigen mode with zero energy, i.e. the 



In the calculation of the other eigen mode around the soliton, we use the soliton form 



obtained in Appendix A. The eigen-mode equation (2.14) is approximated by 



a2 



+ 



+ 



^6 



12 d{Kxy 120 d{Kx)^ 840 d{Kx 



\A(x) 



Ar 



5^ (B.2) 



where we expanded Dj until 6th derivatives with respect to kx. Using (A. 6) with a = 1.70, 
the eigen mode is determined by the numerical diagonalization. There are two additional 
modes with localized character. The eigen values of the localized modes are summarized in 
Table I. They have a good coincidence with the result given by ref.|23]. 
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Figure Captions 



Fig. 1: The loop diagrams which contribute to the effective potential of the lattice 
distortion. These diagrams do not have the momentum transfer between the 
electron and the lattice. 

Fig. 2: The loop diagrams which contribute to the effective potential of the lattice 
distortion. These diagrams have the momentum transfer q between the electron 
and the lattice. 

Fig. 3: The A dependence of the effective potential. The dotted line is the effective 
potential which is given by the Hamiltonian ( 2.14| ). The solid line is the potential 

o 

Fig. 4: The soliton configuration. The dotted line is the solution of the soliton equation 
which is approximated by the second order differential equation ( A.2| ). The solid 



line is the soliton which is given by the variational method. This form is nearly 
equal to the form tanh(x/^ij). 

Fig. 5: The diagrams of the 3 point and 4 point vertices which appear in the Hamilto- 



man 



Fig. 6: The one loop renormalization diagrams for the 3 point vertex and 4 point vertex, 
respectively. 



Fig. 7: The dji ^ dependence of soliton width. The soliton is obtained by the variational 



method using the renormalized Hamiltonian ( 4.27 ), where the trial function is 
chosen to be ( |A.6| ). Since the function ( [A. 61 ) has the form nearly equal to 
tanh(0.567ax/.^o)! the soliton width is given by = ^o/(0-567a). 

Table I: The frequencies of the localized modes around the soliton. The T modes correspond 
to the translational mode of soliton and the B modes come from the third localized mode 
around the soliton. 



Table II : We summarized the values of lo'^ / {2Xu!q'^) for the localized modes of the lattice 
fluctuations around the soliton. They are obtained by the numerical diagonalization of 
( 2.14| ) which is approximated by the differential equation ( |B.2| ). 

Table III : The bare values of the phonon frequencies and the electron-phonon coupling 
constants. These values were used in the analysis of the Raman scattering experiment of 
t-polyacetylene by Vardeny et. alB> 
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